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Previous fits of sterile neutrino dark matter models to cosmological data assumed a peculiar pro- 
duction mechanism, which is not representative of the best-motivated particle physics models given 
current data on neutrino oscillations. These analyses ruled out sterile neutrino masses smaller than 
8-10 keV. Here we focus on sterile neutrinos produced resonantly. We show that their cosmological 
signature can be approximated by that of mixed Cold plus Warm Dark Matter (CWDM). We use 
recent results on ACWDM models to show that for each mass > 2 keV, there exists at least one 
model of sterile neutrino accounting for the totality of dark matter, and consistent with Lyman-a 
and other cosmological data. Resonant production occurs in the framework of the vMSM (the ex- 
tension of the Standard Model with three right-handed neutrinos). The models we checked to be 
allowed correspond to parameter values consistent with neutrino oscillation data, baryogenesis and 
all other dark matter bounds. 



The sterile neutrino is a very interesting Dark Matter 
candidate [1, 2, 3, 4, 5, 6]. The existence of sterile neu- 
trinos (right-handed or gauge singlet) is one of the most 
simple and natural explanations of the observed flavor 
oscillations of active neutrinos. It was observed long ago 
that such particles can be produced in the Early Universe 
through oscillations with active neutrinos [1]. For any 
mass (above ~ 0.4 keV, which is a universal lower bound 
on any fermionic Dark Matter (DM) particle, see [7] and 
references therein) sterile neutrinos produced in this way 
can end up with a correct relic density [1, 2, 4, 8, 9, 10]. 

A single right-handed neutrino would be unable to ex- 
plain the two observed mass splittings between Standard 
Model (SM) neutrinos. Moreover, should this neutrino 
play the role of DM, its mixing with active neutrinos 
would be too small for explaining the observed flavor os- 
cillations [G, 11]. However, in presence of three right- 
handed neutrinos (one for each SM flavor), active neu- 
trino mass splittings and DM may be explained at the 
same time [6]. Moreover, the mass of each sterile neu- 
trino can be chosen below the electroweak scale and ad- 
ditionally explain the matter-antimatter asymmetry of 
the Universe (baryogenesis) [6]. These observations mo- 
tivated a lot of recent efforts for developing this model, 
called the ^MSM [8, 12, 13, 14, 15, 16], and for constrain- 
ing sterile neutrino DM [17, 18, 19]. 

Because of its mixing with flavor neutrinos, this DM 
particle has a small probability of decaying into an ac- 
tive neutrino and a photon of energy E = m s /2 [20], 
producing a monochromatic line in the spectrum of DM 
dominated objects. The corresponding photons flux de- 
pends on the sterile neutrino mass m s and mixing angle 8 



as F ~ 9 2 m 5 s . For each value of the mass and of other pa- 
rameters in the model (sec below) , the angle is fixed by 
the requirement of a correct DM abundance. Combining 
this constraint with bounds on decay lines in astrophys- 
ical spectra allows to put an upper limit on the mass of 
DM sterile neutrinos [3, 5, 17, 18, 19]. 

Within the ^MSM, the relation between m s , 9 and the 
DM abundance can be affected by the presence of a lep- 
ton asymmetry (an excess of leptons over anti-leptons). 
In this case, the production of sterile neutrinos may be 
of the resonant type [2]. The lepton asymmetry re- 
quired for this mechanism to be effective is several or- 
ders of magnitude larger than the baryon asymmetry 
Vb ~ 10~ 10 . In many models of baryogenesis (for a re- 
view see e.g. [21]), both asymmetries are of the same 
order, because they are generated above the electroweak 
scale and sphaleron processes equalize them. Instead, in 
the z/MSM, the lepton asymmetry is generated below the 
electroweak scale, when sphaleron processes are not ac- 
tive anymore [9], As a result, it can be as large as the 
upper limit imposed by Big Bang Nucleosynthesis (BBN) 
and other cosmological constraints (see e.g. [22] and refs. 
therein). Such a large lepton asymmetry is consistent 
with generic values of the parameters of the z/MSM, satis- 
fying current data on neutrino oscillations, cosmological 
requirements (baryogenesis, BBN constraints) and par- 
ticle physics constraints [23]. So, resonant production 
(RP) is a natural way of producing sterile neutrino DM 
in the fMSM. At the same time, most previous constrains 
on sterile neutrino DM assumed non-resonant production 
(NRP) [1, 8]. 

In the NRP case, the comparison of X-ray bounds [17, 
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18, 19] with constraints on DM relic abundance [8] gives 
an upper bound m" RP < 4 keV on the DM sterile neutrino 
mass. In the more effective RP scenario, smaller mixing 
angles are required and the corresponding bound is much 
weaker: m RP < 50 keV [10]. 

The most robust lower bound on the DM mass comes 
from the analysis of the phase space density of compact 
objects, e.g. dwarf spheroidals of the Milky Way halo. 
The universal Gunn-Tremaine bound [24] can be made 
stronger if one assumes a particular primordial phase- 
space distribution function. For NRP sterile neutrinos, 
this leads to m^ RF > 1.8 keV, while for RP particles the 
bound is weaker: m RP > 1 keV [7]. 

An interesting property of sterile neutrino DM with 
keV mass is that it falls in the Warm Dark Matter 
(WDM) category. Lyman-a (Ly-a) forest observations in 
quasar spectra provide strong lower bounds on the mass 
of DM sterile neutrinos produced with the NRP mecha- 
nism [25, 26, 27]. The analysis of SDSS Ly-a data led to 
raf p > 13 keV in Ref. [26], or m™ p > 10 keV in [25]. In 
[27] these bounds were revisited using the same SDSS Ly- 
a data (combined with WMAP5 [28]), but paying special 
attention to the interpretation of statistics in the param- 
eter extraction, and to possible systematic uncertainties. 
It was shown that a conservative (frequentist, 3-cr) lower 
bound is m™ F > 8 keV. The Ly-a method is still un- 
der development, and there is a possibility that some of 
the related physical processes arc not yet fully under- 
stood. However, at this moment, it is difficult to identify 
a source of uncertainty that could give rise to system- 
atic errors affecting the result by more than 30%. Even 
with such an uncertainty, the possibility to have all DM 
in the form of NRP sterile neutrinos is ruled out by the 
comparison of Ly-a results with X-ray upper bounds [27]. 

In the RP case, Ly-a bounds have not been derived 
yet. However, in Ref. [27], a ACWDM model - contain- 
ing a mixture of WDM (in the form e.g. of NRP sterile 
neutrinos) and Cold Dark Matter (CDM) - was analyzed. 
Below we will show that although the phase-space distri- 
bution of RP sterile neutrinos does not coincide exactly 
with such mixed models, some results can be inferred 
from the ACWDM analysis. In particular, we will show 
that for each mass > 2 keV, there is at least one value of 
the lepton asymmetry for which the RP sterile neutrino 
model is fully consistent with Ly-a and other cosmolog- 
ical data (this value of the lepton asymmetry is natural 
within the idVISM). 

Spectra of RP sterile neutrino. DM production in 
the RP scenario occurs in two stages [2, 10]. In pres- 
ence of a lepton asymmetry, the conditions for reso- 
nant oscillations (related to the intersection of disper- 
sion curves for active and sterile neutrinos) are fulfilled 
for temperatures of few hundred MeV. Later, at T ~ 
150 (m^/keV) 1 / 3 MeV, non-resonant production takes 
place. As a result, the primordial velocity distribution of 
sterile neutrinos contains a narrow resonant (cold) com- 
ponent and a non-resonant (warm) one. Its exact form 
can be computed by taking into account the expansion of 
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FIG. 1: Characteristic form of the RP sterile neutrino distri- 
bution function for m" p = 3 keV and various values of the lep- 
ton asymmetry parameter Lq. The spectrum for — 16 (red 
solid line) is shown together with its resonant (dashed) and 
non-resonant (dashed-dotted) components. All these spectra 
have the same shape for q > 3. 



the Universe, the interaction of neutrinos with the con- 
tent of the primeval plasma and the changes in the lepton 
asymmetry resulting from DM production. In this work, 
we used the spectra computed in Ref. [10]. 

Characteristic forms of the spectra are shown in Fig. 1. 
They are expressed as a function of the comoving mo- 
mentum q = pjT v (T v being the temperature of active 
neutrinos) , and depend on the lepton asymmetry param- 
eter [9] L 6 = lQ 6 (n Uc — n,p e )/s (s being the entropy den- 
sity). In rest of this work, the notation M2L25 would 
refer to a model with m RP = 2 keV and Lq = 25. The 
shape of the non-resonant distribution tail depends on 
the mass, but not on Lq. For q > 3, the distribution is 
identical to a rescaled NRP spectrum [8] with the same 
mass (red dashed-dotted line on Fig. 1). We call this 
rescaling coefficient the warm component fraction / NRP . 
For the few examples shown in Fig. 1, the M3L16 models 
corresponds to / NRP ~ 0.12, M3L10 to / NRP ~ 0.53 and 
M3L25 to / NRP ~ 0.60. The maximum of q 2 f(q) for the 
NRP component occurs around q « 1.5 — 2. We define 
the cold component to be the remaining contribution: its 
distribution is given by the difference between the full 
spectrum and the rescaled NRP one (red dashed line on 
Fig. 1), and peaks around (7 rcs ~ 0.25 — 1. Its width and 
height depend on Lq and m RP . 

The DM clustering properties can be characterized 
qualitatively by the particle's free- streaming horizon (see 
e.g. [27] for definition), proportional to its average veloc- 
ity (q)/m. In the RP case, the dependence of the av- 
erage momentum (q) on m RP and Lq is not monotonia 
For a given mass, the RP model departing most from an 
NRP model is the one with the smallest (q) (i.e., with 
the most significant cold component). For each mass, 
there is indeed a value of Lq minimizing (q), such that 
(q}min ~ 0.3(g) NRP (cf. [10]). This minimum corresponds 
to lepton asymmetries that are likely to be generated 
within the z/MSM. For instances, for m RP = 2, 3 or 
4 keV, the spectra with the smallest average momentum 
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are M2L25, M3L16 and M4L12, all having / NRP < 0.2. 

For a quantitative analysis, we computed the power 
spectrum of matter density perturbations P KP (k,z) for 
these models. The standard software (i.e. CAMB [29]) is 
not immediately appropriate for this purpose, as it only 
treats massive neutrinos with a Fermi-Dirac primordial 
distribution. To adapt it to the problem at hand, we 
modified CAMB so that it could take arbitrary spectra as 
input data files. We analyzed the spacing in momentum 
space needed in order to obtain precise enough results, 
and implemented explicit computations of distribution 
momenta in CAMB. We cross-checked our results by mod- 
ifying another linear Boltzmann solver - CMBFAST [30], 
implementing a treatment of massive neutrinos with ar- 
bitrary analytic distribution function. 

To separate the influence of primordial velocities on the 
evolution of density perturbations from that of cosmolog- 
ical parameters, it is convenient to introduce the transfer 

function (TF) T(k) = [P HP (fc)/P A CDM(fc)] 1/2 ■ Figure 2 
shows the transfer function of the models M3L16 and 
M4L12. The TF becomes smaller than one above the 
wave number associated with the free-streaming horizon 
today, fc FSH w 0.5 (mf /l keV) hMpc' 1 (c.f. [27]). We 
see that for a large range of k values above fc FSH , roughly 
k < 5/c F sh, the transfer function T RP (fc) is very close to 
Tacwdm(^) for the same mass and warm component frac- 
tion F WDM = / NRP . On smaller scales, T RP (k) decreases 
faster, since the cold component of RP sterile neutrinos 
also has a non-negligible free-streaming scale. For all 
values of the mass studied in this work, mf > 2 keV, 
the discrepancy appears above 5 ft,/Mpc (vertical line 
in Fig. 2), i.e. above the maximum scale in the three- 
dimensional power spectrum to which current Ly-a data 
are sensitive. Hence, for the purpose of constraining RP 
sterile neutrinos with Ly-a data, it is possible to use 
the results obtained in the ACWDM case. In Ref. [27], 
we presented the results of a WMAP5 plus SDSS Ly-a 
data analysis for ACWDM models with mf > 5 keV. 
In Fig. 3, we show the Bayesian credible region for the 
mass and the warm component fraction, now extended 
till mf p = 2 keV [34]. 

Fig. 2 demonstrates that for the models M3L16 and 
M4L12, the function T RP (k) lies above TacwdmCO for 
the same masses and F WDM = 0.2 (at least, in the range 
of wave numbers probed by Ly-a data): so, it can only be 
in better agreement with cosmological data. We checked 
that the same is true for M2L25. However, ACWDM 
models with mf = 2, 3, 4keV and F WDM = 0.2 are within 
the 2-er contour of Fig. 3. We conclude that M2L25, 
M3L16 and M4L12 are clearly allowed by the data. For 
larger mass (and still minimal (q)), the free-streaming 
horizon is smaller, and agreement with observations can 
only become easier. Therefore, we see that for each mass 
?n RP > 2 keV there exists at least one value of the lepton 
asymmetry for which RP sterile neutrinos are perfectly 
compatible with WMAP5 and SDSS Ly-a data. 

Fig. 4 shows the range of masses and mixing angles 
consistent with constraints from phase-space density ["i 
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FIG. 2: TFs for the models M4L12 (top) and M3L16 (bot- 
tom), together with CWDM spectra for the same mass and 
Fwdm ^ 0.15 or 0.2. 
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FIG. 3: 1 and 2-a bounds from WMAP5 and SDSS Ly-a data 
on ACWDM parameters. Red points correspond to approx- 
imations for the models M4L12, M3L16, M2L25. Results to 
the left of black vertical lines were already reported in [27]. 



(left shaded region), from X-rays [17, 19] (upper right 
corner, shaded in red) and providing the correct DM 
abundance (curves between the lines "NRP" and L™ ax : 
from top to bottom L 6 = 8, 12, 16, 25, 70, 250). The black 
dashed line shows approximately the RP models with 
minimal (q) for each mass, i.e., the family of models with 
the largest cold component. We have seen that all black 
filled circles along this line and with mf > 2 keV are 
compatible with Ly-a bounds. In addition, those with 
mf < 4 keV arc also compatible with X-ray bounds 
(this conclusion does not change with the new results of 
Ref. [31]). Note that above 4keV, Ly-a data allows in- 
creasingly high WDM fractions, so that agreement with 
both Ly-a and X-ray bounds can be maintained with 
larger values of L§. This is very clear e.g. for the models 
M10L25, M10L16 and M10L12, allowed by X-ray data 
(open circles on Fig. 4), and consistent with Ly-a data 
since for mf = 10 keV, up to 100% of WDM is allowed 
at the 2-cr level (c.f. Fig.3). 

In conclusion, we showed in this work that sterile neu- 
trino DM with mass > 2keV is consistent with all existing 
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FIG. 4: Region of masses and mixing angles for RP sterile 
neutrinos consistent with existing constraints. 

constraints. A sterile neutrino with mass ~ 2 keV is an 
interesting WDM candidate, as it may affect structure 
formation on galactic scales. This range of masses and 
corresponding mixing angles is important for laboratory 
and astrophysical searches. 



To determine the precise shape of the allowed param- 
eter range (which may continue below 2 keV, see Fig. 4), 
one should perform specific hydrodynamical simulations 
in order to compute the flux power spectrum on a grid 
of (m" Lq) values, and compare with Ly-a data. We 
leave this for future work. 

The z/MSM does not require new particles apart from 
the three sterile neutrinos. Extensions of this model may 
include a scalar field providing Majorana masses to ster- 
ile neutrinos via Yukawa couplings [32, 33]. Then, sterile 
neutrino DM can also be produced by the decay of this 
scalar field, and also contain a cold and a warm compo- 
nent. We expect a similar range of masses and mixing 
angles to be allowed by Ly-a data. The quantitative 
analysis of this model is also left for future work. 
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Thermal velocities of NRP neutrinos may become im- 
portant below 5 keV [27]. To be conservative, this mass 
region was excluded from the results of [27], but we put it 
back here. We see that 1 and 2a contours become nearly 
horizontal for m < 5 keV (c.f. Fig. 3). We do not expect 
thermal velocities to affect this conclusion. For detailed 
bounds on RP sterile neutrinos, we plan to explore the 
precise dependence of N-body simulation results on ther- 
mal velocities elsewhere. 



